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ABSTRACT 

Sub-millimetre observations of the William Herschel Deep Field (WHDF) using 
the Large Apex Bolometer Camera (LABOCA) revealed possible sub-mm counterparts 
for two X-ray absorbed quasars. The primary aim here is to exploit Expanded Very 
Large Array (EVLA) radio continuum imaging at 8.4 GHz to establish the absorbed 
quasars as radio/sub-mm sources. The main challenge in reducing the WHDF EVLA 
data was the presence of a strong 4C source at the field edge. A new calibration 
algorithm was applied to the data to model and subtract this source. The resulting 
thermal noise limited radio map covers a sky area which includes the 16' x 16' Extended 
WHDF. It contains 41 radio sources above the 4tr detection threshold, 17 of which 
have primary beam corrected flux densities. The radio observations show that the 
two absorbed AGN with LABOCA detections are also coincident with radio sources, 
confirming the tendency for X-ray absorbed AGN to be sub-mm bright. These two 
sources also show strong ultraviolet excess (UVX) which suggest the nuclear sightlinc 
is gas- but not dust-absorbed. Of the three remaining LABOCA sources within the 
~ 5' half-power diameter of the EVLA primary beam, one is identified with a faint 
nuclear X-ray /radio source in a nearby galaxy, one with a faint radio source and one 
is unidentified in any other band. 

More generally, differential radio source counts calculated from the beam-corrected 
data are in good agreement with previous observations, showing at S < 50/iJy a sig- 
nificant excess over a pure AGN model. In the full area, of ten sources fainter than 
this limit, six have optical counterparts of which three are UVX (i.e. likely quasars) 
including the two absorbed quasar LABOCA sources. The other faint radio counter- 
parts are not UVX but are only slightly less blue and likely to be star-forming/merging 
galaxies, predominantly at lower luminosities and redshifts. The four faint, optically 
unidentified radio sources may be either dust obscured quasars or galaxies. These high 
redshift obscured AGN and lower redshift star-forming populations are thus the main 
candidates to explain the observed excess in the faint source counts and hence also the 
excess radio background found previously by the Absolute Radiometer for Cosmology, 
Astrophysics and Diffuse Emission (ARCADE2) experiment. 

Key words: galaxies: high-redshift - quasars: general - radio continuum: galaxies - 
techniques: interferometric 
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1 INTRODUCTION 

The William Herschel Deep Field (WHDF; 7x7 arcminutes 
at OOh 22m 30s +00° 21m 00s) has been comprehensively 
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targeted by many authors using many observatories for over 
two decades (Metcalfe et al., 1991, 1995, 2001, 2006; Mc- 
Cracken et al., 2000). It is one of the deepest survey fields in 
terms of optical and near infrared photometry, having more 
than 150 hours of CCD imaging on the 4.2 m William Her- 
schel Telescope (WHT) and 3.9 m telescope at Kitt Peak Na- 
tional Observatory, resulting in magnitude limits of U = 27, 
B = 28.2, R = 27 and I = 25.5. The magnitude apertures 
vary slightly with band / seeing but are « l."3 radius for the 
faintest stars. All magnitudes quoted are in the Vega system. 
More details are given by Metcalfe et al. (2001). There is also 
deep H band imaging following 30 hours of exposure with the 
Calar Alto 3.5 m telescope (Vallbe-Mumbru, 2004). Deep, 
high-resolution I-band imaging is available through Hub- 
ble Space Telescope (HST) Advanced Camera for Surveys 
(ACS) imaging (Bohm & Zeigler, 2006). Observations with 
the Low Dispersion Survey Spectrograph 2 (LDSS-2) spec- 
trograph on the 6.5 m Magellan telescope (Vallbe-Mumbru, 
2004) and Very Large Telescope (VLT, 8.2 m) observations 
using FORS2 (Bohm & Zeigler, 2006) have yielded ~200 
spectroscopic redshifts in the field. 

Moving away from optical and infrared wavelengths, the 
WHDF has deep X-ray imaging via 75 ks of Chandra ACIS- 
I data (Vallbe-Mumbru, 2004), and Bielby et al. (2012) 
recently obtained sub-mm observations of the field using 
the Large Apex Bolometer Camera (LABOCA), reaching 
a depth of 1 mjy at a wavelength of 870 um and resulting 
in the detection of eleven sources. 

This paper presents the latest addition to the multi- 
wavelength coverage of the WHDF, namely 8.4 GHz ob- 
servations using the Expanded Very Large Array (here- 
after EVLA, since renamed the Karl G. Jansky Very Large 
Array). These observations were motivated by the recent 
LABOCA data with the high resolution afforded by the ra- 
dio observations being required to confirm whether any of 
the bright sub-mm sources in the field were associated with 
the significant sample of both absorbed and unabsorbed 
quasars (QSOs) also present in the WHDF. In particular, 
Bielby et al. (2012) suggested two X-ray absorbed high red- 
shift active galactic nuclei (AGN), one of which is a Type 2 
quasar showing only narrow lines in the optical, may be asso- 
ciated with LABOCA sources, in line with previous sugges- 
tions of a connection between X-ray absorbed AGN and sub- 
mm emission (Gunn 1999; Page et al. 2004; Hill & Shanks 
2011). It is therefore important to establish the reality of the 
association between the sub-mm and these absorbed AGN 
as far as possible, using these radio data. 

The WHDF is hitherto largely unexplored at radio 
wavelengths due to the presence of a strong 4C source 6 
arcminutes south of the field centre, however this has been 
cleanly subtracted from the observations described in this 
paper using a new calibration algorithm. A description of 
the observations and the calibration process follows in Sec- 
tion [2] The radio map and source catalogue are presented 
in Section [3] Section U matches the radio sources with coun- 
terparts at other wavebands. These results are discussed in 
Section [5j and concluding remarks are made in Section [6] 



2 OBSERVATIONS AND DATA REDUCTION 

The WHDF was observed with the EVLA in its most com- 
pact D-configuration between 28 March 2010 and 15 June 
2010. A total of 35 hours were awarded under the Open 
Shared Risk Observing program, of which 30 hours were ob- 
served. The observations were recorded in 25 independent 
scheduling blocks (SBs) with durations of either 1.5 or 3.5 
hours, resulting in 25 sets of visibility data. 

The X-band receivers were used and the correlator was 
configured to deliver two spectral windows centred at 8.396 
and 8.524 GHz, each providing 64 x 2 MHz channels giving 
256 MHz of contiguous frequency coverage. For this contin- 
uum experiment an 8:1 frequency averaging was applied to 
make the data volume more manageable. No time averag- 
ing was applied, so the standard 1 second integration time 
was preserved. Each of the 25 Measurement Sets were then 
split further in order to separate each of the two spectral 
windows, resulting in a total of 50 unique data sets to be 
independently edited and calibrated. 

Editing of the data was performed using the PlotMS 
tool within the NRAO CASAp package. In general the data 
were free from corruptions save for the occasional bursts 
of radio frequency interference which were simple to locate 
and excise. Two SBs exhibited severe amplitude drifts on a 
majority of baselines and in both spectral windows and were 
discarded outright. 

For flux calibration a single observation of J0137+3309 
(3C48) was performed at the start of each SB, and the phase 
calibrator J0022+0014 (4C 00+02; 0.64 Jy at X-band) was 
observed for 40 seconds for every 9 minutes of target obser- 
vation. 

Standard calibration procedures were followed using 
CASA. The flux scale for the amplitude calibrator was set 
against a model image as 3C48 is slightly resolved in X-band 
observations, and this source was also used to calibrate the 
bandpass. The per-antenna complex gain solutions were gen- 
erated for the phase calibrator source and then interpolated 
across the target field. A Python script was constructed to 
apply this standard procedure to each of the 50 flagged Mea- 
surement Sets automatically. Plots of the gain solutions were 
generated and inspected in order to ensure successful cali- 
bration. Phase stability was excellent throughout. As an ad- 
ditional diagnostic, calibrated maps of the phase calibrator 
and target field were also generated by the script. 

Inspection of the images of the target field immediately 
shows why the WHDF is a challenge for radio interferome- 
ter observations, as the phase calibrator is situated approxi- 
mately 6 arcminutes south of the target field and completely 
dominates the radio emission in each map. Furthermore, as 
with many well-studied multiwavelength fields, the WHDF 
is situated close to the celestial equator (Declination +20 
arcminutes). When observing the celestial equator with an 
interferometer the locus traced by a given baseline in the uv 
plane as the Earth rotates exhibits a strong east- west bias. 
The Fourier transform of the uv coverage determines the 
point-spread function (PSF, also known as the dirty beam) 
of the observation, thus a uv plane sampling pattern that 
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is dominated by east-west structure results in strong north- 
south sidelobes in the PSF. As can be seen in the upper 
panel of Figure [1] the sidelobes associated with the phase 
calibrator completely cover the target field, and these ef- 
fects must be mitigated in order to obtain a scientifically 
useful radio map. 

2.1 Subtracting the phase calibrator from the 
target field 

Straightforward multifrequency synthesis deconvolution of 
the phase calibrator from the target field using the CLEAN 
algorithm did not yield the result that one might hope for, 
and residual sidelobe structure still swamped the target 
field: the upper panel of Figure [T] actually shows the cali- 
brated target field after deconvolution has been attempted. 
Similarly, subtraction of a visibility model derived from the 
clean components followed by imaging of the residual data 
set resulted in a map that was still corrupted by strong 
north-south emission associated with the phase calibrator. 

Examining a model of the EVLA primary beam at 8.4 
GHz (Walter Brisken, private communication; Figure IA2|) 
reveals that when the array is pointing at the target field 
the phase calibrator lies at the boundary of the first null and 
the first sidelobe, with its radial separation from the point- 
ing centre changing with frequency. The first sidelobe of the 
EVLA beam has four-way azimuthal structure caused by the 
Cassegrain optics of the dishes, and the azimuth-elevation 
mount of the EVLA dishes causes this pattern to rotate on 
the sky as a source is tracked. Since the primary beam can 
to first order be thought of as the gain of the instrument as a 
function of direction, these effects conspire to impart appar- 
ent temporal variability to the confusing source (as well as 
significant spectral effects). Sources which vary significantly 
on timescales shorter than the observation have corrupted 
PSFs associated with them, and this is why traditional de- 
convolution and model subtraction cannot cope. 

A calibration scheme which can solve for these 
direction-dependent effects must be employed, and for these 
observations the differential gains algorithrrQ (Smirnov, 
2011b) was used, implemented using the Calico frame- 
work within the MeqTree^f] software package (Noordam & 
Smirnov, 2011). The algorithm works by solving for addi- 
tional complex gain terms against an assumed sky model 
for a subset of sources. In the case of these observations the 
solutions are generated for a sky model consisting solely of 
a point source at the location of the phase calibrator. Since 
this source dominates the radio emission, the gain solutions 
will also be dominated by corruptions in the direction of that 
source. There is no need to provide a stringent measure of 
the source flux for the sky model, as any errors in this pa- 
rameter will be subsumed into the gain solutions, however 
the flux in the sky model was fixed at 5 mjy, as measured 
from the image generated following the initial calibration 
pass. Note the high attenuation of this source by comparing 
this measured value to its intrinsic flux measured by imaging 
the phase calibration scans. 

2 Also known as "The Flyswatter" for reasons that will become 
obvious. 

3 http://www.astron.nl/meqwiki 



Data are averaged over five minute intervals and in 
pairs of frequency channels (i.e. four blocks across the aver- 
aged band). A solution is generated for each of these time- 
frequency tiles. This has the effect of reducing the degrees of 
freedom in the fit, whilst simultaneously time-smearing out 
contributions from the rest of the field as well as accounting 
for any spectral behaviour in either the direction-dependent 
gain corruptions or the source itself. In this regime, the use of 
the differential gains algorithm is analogous to the more gen- 
erally employed "peeling" algorithm (e.g. Noordam, 2004), 
however the advantage of the differential gains technique for 
more general applications is one of flexibility: it can generate 
additional gain terms for many individual sources simulta- 
neously, whilst simultaneously solving on shorter timescales 
for the traditional complex receiver gains derived from an 
all-inclusive sky model. Peeling generally requires a cum- 
bersome iterative approach whereby dominating sources are 
treated in order of brightness. 

The best fit visibility model derived from this process 
was subtracted from the observed visibilities, and the resid- 
ual data were imaged. Following this procedure the confus- 
ing source is completely removed leaving no residual emis- 
sion above the noise. As a quantitative example of the suc- 
cess of this of this technique, the root-mean-square (rms) 
background level in a map formed from one of the spectral 
windows of one of the 3.5-hour scheduling blocks was 21 /ijy 
/ beam (as measured away from the dominating residual 
north-south structure) following either deconvolution of the 
confusing source, or subtraction of an inverted clean com- 
ponent model from the visibilities. This value dropped to 
12 /ijy / beam following the subtraction of the confusing 
source using the differential gains algorithm. The radio im- 
ages corresponding to these "before" and "after" scenarios 
can be seen in Figure [T] with further details provided in the 
caption. 

Following the discarding of the two corrupted SBs, the 
first step in the calibration produced 46 Measurement Sets 
which were then split in order to contain only calibrated 
visibilities from the target field. MeqTrees can be operated 
non-interactively, allowing the source subtraction process to 
be scripted and applied automatically to each of these Mea- 
surement Sets. Again, calibrated maps were produced for 
each individual Measurement Set and the differential gain 
solutions were examined to check for problems. The gain so- 
lutions themselves encode information as to the origin of the 
direction-dependent corruptions, and these are discussed in 
the Appendix. 



2.2 Concatenation and imaging 

The calibrated data with the confusing source removed were 
concatenated into a final Measurement Set which was im- 
aged to form the final map. At this stage, only 23 of the 
46 viable Measurement Sets were included. The ones which 
were not included were shorter duration runs at low eleva- 
tions. The observations had relaxed hour angle constraints 
when they were scheduled, and the inclusion of these data 
did not improve the quality or depth of the final map, which 
is presented in Section \3. II 
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Figure 1. Radio images of the target field generated from a sin- 
gle 3.5-hour Measurement Set. The upper panel shows the image 
that results from the initial complex gain calibration and a multi- 
frequency synthesis deconvolution of the phase calibrator / con- 
fusing source using the CLEAN algorithm. The time-dependent 
gains imparted by the primary beam of the array in that direction 
corrupt the PSF associated with this source and even following 
deconvolution the residual sidelobe structure completely covers 
the science targets. The lower panel shows the map that results 
following subtraction of the confusing source using the differential 
gains algorithm (Smirnov, 2011b). The position of the confusing 
source is marked with the cross. Radio sources in the target field 
are now visible. The contour levels in these maps are the same, 
and begin at 3cr (where cr = 12 ^tJy, the background rms of the 
lower map) and increase in multiples of 2. A single, dashed neg- 
ative countour is included with a value of -3<r. 



3 RESULTS 

The following two sections present the final radio map of the 
WHDF at 8.4 GHz and the source catalogue derived from 
it. 

3.1 Radio map 

Figure [2] shows the final, deconvolved 8.4 GHz radio image 
covering the WHDF. Deconvolution was performed using 
the CASA clean task in multi-frequency synthesis mode with 
w-term correction. The restoring beam is a circular Gaussian 
with a full- width at half-maximum spanning 8 arcseconds, 
as shown by the filled circle in the lower left-hand corner of 
the map. The rms of the background noise level in this map 
is 2.5 /xJy. This is consistent with the theoretical thermal 
noise level expected in this observation to within 3%. The 
contours on Figure [2] begin at 3a and increase in multiples 
of V2. There is a single dashed negative contour at -3a. 

The area of sky covered by this image is 22.2 x 22.2 
arcminutes, although most of the multiwavelength observa- 
tions occupy a 7 x 7 arcminute patch near the map centre. 
This is well-matched to the primary beam of the EVLA an- 
tennas, the approximate half-power point of which is shown 
by the large dashed circle in Figure [2] Imaging over this ex- 
tended area in the radio was done to search for bright radio 
sources in the 16 x 16 arcminute area covered by the shal- 
lower Extended WHDF observations. The position of the 
phase calibrator source (4C 00+02) that was confusing the 
central region before subtraction is also marked. No resid- 
ual emission is present above the noise. Also noteworthy 
is the spatially extended radio source towards the western 
edge of the map (WHDF-EVLA-3 in Table [TJ. This source 
is present in the Sloan Digital Sky Survey (SDSS; Abazajian 
et al., 2009) which lists it as a galaxy at z = 0.2609. 

The crosses on Figure [2] show the positions of the 
spectroscopically-confirmed quasars in the WHDF (Vallbe- 
Mumbru, 2004). The small circles show the locations of the 
LABOCA-detected sub-mm sources in the field, with the 
size of the circles being the positional uncertainties in the 
detections. The nomenclature of these two classes of sources 
on Figure [2] and in the discussion that follows is a concise 
version of that employed by Bielby et al. (2012): three digit 
identifications are confirmed quasars and two digit identifi- 
cations correspond to LABOCA sources (prefixed by WHD- 
FCH and WHDF-LAB respectively in the aforementioned 
article). Radio sources are referred to using the full WHDF- 
EVLA prefix throughout, with the exception of the labels 
on Figure [3] This covers the same sky area as Figure [2] but 
shows the location of each source according to its identifi- 
cation as listed in Table [T] Figure [3] also shows the central 
and extended regions of the WHDF. Higher magnification 
thumbnails of individual radio sources can be found in Fig- 
ures[4]and[5] along with images of the counterparts identified 
at other wavebands (see Section I4.2[l . 

3.2 Source catalogue 

The radio map presented in Section 13.11 was searched for 
peaks of emission exceeding 4a (= 10 pJy) using the SAD 
task within the Astronomical Image Processing System 
(AIPS; Greisen, 2003) package. The resulting catalogue was 
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Figure 2. Radio contour map of the William Herschel Deep Field at 8.4 GHz. The rms noise (<r) in the map is 2.5 fi3y / beam. 
Contours begin at 3cr and increase in multiples of x/2. There is a single dashed negative contour at -3<r. The numbered crosses show 
the spectroscopically confirmed quasars in the WHDF (Vallbe-Mumbru, 2004) and the numbered circles indicate the locations of the 
LABOCA-detected sub-mm sources (Bielby et al., 2012). The numbers associated with these two classes of sources are consistent with 
those used by Bielby et al. (2012), see Section l3.1l for details. The size of the circles used to show the locations of the LABOCA sources 
indicates the positional uncertainty. The large dashed circle shows the approximate half-power point of the EVLA primary beam. The 
position of the phase calibrator 4C 00+02 is also marked although no detectable emission remains. Note also the extended SDSS source 
to the west. The filled circle in the lower left-hand corner of the plot shows the 8 arcsecond extent of the circular restoring beam applied 
following deconvolution. 



manually pruned as two of the sources (WHDF-EVLA-3 and 
WHDF-EVLA-6) were resolved into multiple Gaussian com- 
ponents. In these cases the flux was determined by summing 
up the components. The 22.2 x 22.2 arcminute map shown 
in Figure [2] contains 41 discrete sources with flux densities 
exceeding 4a and the properties of these are listed in de- 
scending order of brightness in Table [T] 



Of the sources detected in the map, 17 of them were 
within the region where the gain of the main lobe of the 
EVLA primary beam exceeded 0.2. These are indicated by 
the circles in Table [T] (see also Figure [3]). An analytic expres- 
sion is available to correct for the beam gain of the (E)VLA 
antennas as per the AIPS task PBCOR, and this has been 
applied to the 17 sources in the main lobe. Thus bold-face 
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type in Table Q] represents intrinsic flux values and plain 
type entries have apparent flux values. The measured sep- 
aration from the pointing centre that was used to derive 
the primary beam correction factor for each source is also 
listed. Uncertainties in the intrinsic fluxes include contribu- 
tions from errors in the fit to the main lobe of the beam. 

The subtracted phase calibrator (0.64 Jy intrinsic flux) 
is not included in this table. If one were interested in this 
source it would be more useful to simply image the phase cal- 
ibrator scans during which the array was pointing directly at 
it rather than trying to correct for its complicated behaviour 
in the target scans. 



4 SUB-MM AND OPTICAL/NEAR- INFRARED 
RADIO SOURCE COUNTERPARTS 

4.1 Sub-mm counterparts 

High redshift extragalactic sources which are luminous at 
sub-mm wavelengths are generally thought to be driven by 
high star formation rates, with the high sub-mm flux in 
the system arising due to the reprocessing of the intense 
stellar radiation by large amounts of dust. The low tem- 
perature of the dust is thought to favour such a heating 
source rather than AGN. However the high redshift sub-mm 
sources clearly have the bolometric luminosities of AGN. 
Moreover, explaining the observed sub-mm source counts as 
being of purely star-formation origin in a Cold Dark Mat- 
ter cosmology conte xt requires invocat ion of top heavy ini- 
tial mass functions (Ba ugh et all 12005). Recently, evidence 
has mounted for AGN ma king significant contributions to 
sub-m m source counts. fe.g. lLutz et aljfeoid : IHiII fc Shanks! 
2011a) following earlier identifications of AGN in sub-mm 
sources (eg Ivison et al. 1998, Brandt et al. 2001) 

These findings motivated the 870 /jm LABOCA ob- 
servations of Bielby et al. (2012) as the WHDF conve- 
niently contains 15 spectroscopically confirmed quasars, four 
of which are classified as obscured via the assumption that 
their high X-ray hardness ratirjf] is caused by the presence 
of large amounts of absorbing hydrogen. Eleven sub-mm 
sources were detected in the LABOCA observations, two 
of which appeared to be associated with heavily obscured 
quasars. None of the unobscured quasars were near to robust 
sub-mm counterparts, although a stacking analysis revealed 
that absorbed X-ray quasars showed the most significant 
sub-mm emission. 

Since sub-mm observations of this kind have relatively 
low resolution, radio observations have been used to de- 
tect counterparts for sub-mm sources by several authors 
(e.g. Ivison et al., 2002), with the improved positional ac- 
curacy afforded by the radio observations then being used 
to confirm counterparts at other wavebands. Recalling the 
naming scheme defined in Section 13.11 the EVLA obser- 
vations detect radio counterparts at the positions of four 
sub-mm sources: 06, 11, 05 and 02 corresponding to radio 
sources WHDF-E VL A- 1 , WHDF-EVLA-6, WHDF-EVLA-8 
and WHDF-EVLA-32. Radio sources WHDF-EVLA-6 and 

4 If H and S represent photon counts in the hard (2-8 keV) and 
soft (0.5-2 keV) X-ray bands then the hardness ratio is defined 
as (H - S) I (H + S). 



WHDF-EVLA-8 correspond to the X-ray absorbed quasars 
007 {z = 1.33) and 008 (z = 2.12) respectively. 

The reliability of the radio and sub-mm sources be- 
ing counterparts (and not chance alignments) is evaluated 
he re using the c o rrecte d Poisson probability (P) as used 
bv lDownes et all (|l986t ). Based on the number densities of 
radio sources at the fluxes of those presented here, the prob- 
abilities that radio sources WHDF-EVL A- 1 , WHDF-EVLA- 
6, WHDF-EVLA-8, and WHDF-EVLA-32 are chance align- 
ments with the LABOCA sub-mm sources are P = 0.003, 
P = 0.023, P = 0.006 and P = 0.032 respectively. These 
probabilities of < 3% indicate that the radio and sub-mm 
signals do indeed originate from the same sources. 

Summarising, of the 5 sub-mm sources in the central 
WHDF area surveyed by EVLA, 4 are radio sources and 
3 were already known to be candidate X-ray sources. The 
EVLA observations confirm the positional coincidence of 
the three sub-mm source with X-ray sources on the as- 
sumption that the LABOCA and EVLA sources are the 
same. One further LABOCA source is close to an EVLA 
source (02/ WHDF-EVLA-32) but neither source is close to 
an optical/near-infrared object. The final sub-mm source 
(04) is unidentified in either radio or optical/near-infrared. 
Thus three out of five sub-mm sources are AGN as indi- 
cated by their X-ray properties and the other two are blank 
in the optical/near-infrared. Two of the three X-ray sources 
are hard X-ray sources and likely to be cold gas absorbed. 
Although the numbers of sub-mm sources are small, the frac- 
tion that are X-ray sources is high and so is the gas-absorbed 

frac tion. 

Hill fc Shanks! (|2011al ). following iGunn fc Shanks! 

(1999), suggested that the FIR flux re-radiated by the dust 
might be proportional to the amount of X-ray radiation 
absorbed and in this case the obscured AG N contribution 
to the s ub-mm background mig ht reach 40%. lfflll fc Shanks! 
l|2011al ) found strong evidence for this through statistical 
analyses of the Extended Chandra Deep Field South 
(ECDFS) X-ray and sub-mm data. Aided by the EVLA 
confirmations of the identification of the 2 absorbed X- 
ray quasars with the sub-mm sources means that there 
is now further evidence for this hypothesis in the new 
WHDF/EVLA dataset. These points will be revisited in 
the discussion in Section [5] 

4.2 Optical/near-infrared counterparts and 
Individual radio sources 

This section presents possible (near-)infrared, visible and 
ultraviolet counterparts for the radio sources and discusses 
noteworthy individual sources. Cut-out images from exist- 
ing observations are plotted for each radio source in the 
central 7' x 7' area of the WHDF in Figure |4j and in the 
extended area in Figure [5] Photometric bands U (365 nm), 
B (445 nm), R (658 nm), I (806 nm) and H (1630 nm) are 
shown, with the telescope used indicated above each col- 
umn. The final column shows the radio image with over- 
laid contours. The base contour level is 2<r and increases in 
multiples of \f2. Each thumbnail spans 25 arcseconds. The 
nearest counterparts to the radio sources are listed in Ta- 
ble [2] Photometric colours are listed, as are redshifts, where 
available. 

WHDF-EVLA-1: The brightest radio source (aside from 
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ID Right Ascension Declination 

(J2000) (J2000) 






WHDF 


EVLA- 


1 


Oh 


22m 


31 


82s 


+0d 


21m 


28.03s 





WHDF 


EVLA- 


2 


Oh 


22m 


15 


86s 


+0d 


21m 


28.32s 




WHDF 


EVLA- 


3 


Oh 


21m 


50 


09s 


+0d 


22m 


3.15s 





WHDF 


EVLA- 


4 


Oh 


22m 


31 


07s 


+0d 


22m 


24.81s 





WHDF 


EVLA- 


5 


Oh 


22m 


29 


15s 


+0d 


18m 


19.78s 





WHDF 


EVLA- 


6 


Oh 


22m 


24 


89s 


+0d 


20m 


10.21s 




WHDF 


EVLA- 


7 


Oh 


22m 


39 


97s 


+0d 


18m 


38.2s 





WHDF 


EVLA- 


8 


Oh 


22m 


22 


85s 


+0d 


20m 


13.98s 





WHDF 


EVLA- 


9 


Oh 


22m 


32 


6s 


+0d 


23m 


28.48s 





WHDF 


EVLA- 


10 


Oh 


22m 


20 


53s 


+0d 


19m 


1.52s 




WHDF 


EVLA- 


11 


Oh 


21m 


55 


42s 


+0d 


25m 


41.58s 




WHDF 


EVLA- 


12 


Oh 


22m 


10 


95s 


±0d 


19m 


50.99s 




WHDF 


EVLA- 


13 


Oh 


22m 


54 


37s 


+0d 


30m 


39.98s 




WHDF 


EVLA- 


11 


Oh 


22m 


21 


83s 


+0d 


16m 


47.18s 





WHDF 


EVLA- 


15 


Oh 


22m 


30 


09s 


+0d 


19m 


41.63s 





WHDF 


EVLA- 


16 


Oh 


22m 


14 


93s 


+0d 


21m 


4.12s 




WHDF 


EVLA- 


17 


Oh 


22m 


42 


96s 


+0d 


19m 


56.28s 




WHDF 


EVLA- 


18 


Oh 


22m 


37 


79s 


+0d 


23m 


36.78s 





WHDF 


EVLA- 


19 


Oh 


22m 


22 


66s 


+0d 


20m 


42.82s 




WHDF 


EVLA- 


20 


Oh 


22m 


53 


44s 


+0d 


16m 


57.03s 




WHDF 


EVLA- 


21 


Oh 


22m 


23 


89s 


+0d 


11m 


14.61s 




WHDF 


EVLA- 


22 


Oh 


21m 


58 


45s 


+0d 


22m 


41.49s 





WHDF 


EVLA- 


23 


Oh 


22m 


32 


02s 


+0d 


18m 


11.78s 




WHDF 


EVLA- 


24 


Oh 


21m 


47 


32s 


+0d 


15m 


59.91s 





WHDF 


EVLA- 


25 


Oh 


22m 


36 


86s 


+0d 


20m 


55.57s 




WHDF 


EVLA- 


26 


Oh 


22m 


13 


91s 


+0d 


18m 


5.77s 





WHDF 


EVLA- 


27 


Oh 


22m 


16 


16s 


+0d 


19m 


46.78s 




WHDF 


EVLA- 


28 


Oh 


22m 


15 


75s 


+0d 


31m 


11.39s 





WHDF 


EVLA- 


29 


Oh 


22m 


21 


64s 


+0d 


21m 


49.48s 




WHDF 


EVLA- 


30 


Oh 


22m 


31 


64s 


+0d 


31m 


25.44s 




WHDF 


EVLA- 


31 


Oh 


22m 


59 


13s 


+0d 


11m 


20.34s 





WHDF 


EVLA- 


32 


Oh 


22m 


28 


25s 


+0d 


21m 


47.24s 




WHDF 


EVLA- 


33 


Oh 


21m 


43 


99s 


±0d 


18m 


5.18s 




WHDF 


EVLA- 


34 


Oh 


21m 


43 


24s 


+0d 


11m 


24.65s 




WHDF 


EVLA- 


35 


Oh 


22m 


59 


7s 


+0d 


12m 


34.81s 




WHDF 


EVLA- 


36 


Oh 


21m 


59 


92s 


+0d 


20m 


41.05s 




WHDF 


EVLA- 


37 


Oh 


21m 


47 


34s 


+0d 


12m 


23.71s 




WHDF 


EVLA- 


38 


Oh 


22m 


30 


07s 


+0d 


30m 


2.58s 





WHDF 


EVLA- 


39 


Oh 


22m 


30 


84s 


+0d 


22m 


53.9s 




WHDF 


EVLA- 


10 


Oh 


22m 


48 


38s 


+0d 


14m 


49.05s 




WHDF 


EVLA- 


11 


Oh 


22m 


11 


57s 


+0d 


25m 


41.96s 





J_ 1 LLJ\- Ll^llol L y 


R p rl 1 11 1; 

± V cXLl 1 U.O 


Beam-corrected 


ell LoLt_. ; 


(11 Iv") 


I ell 1 


flux density (/^Jy) 


0.1,0.1) 


164.42 (±6.0) 


1.389 


201.23 (±7.55) 


0.08,0.09) 


106.23 (±4.0) 


2.884 


275.71 (±18.4) 


0.64,0.77) 


87.58 (±10.0) 


9.33 


N/A 




0.13,0.12) 


70.88 (±4.0) 


1.939 


106.16 (±6.84) 


0.24,0.26) 


43.85 (±5.0) 


2.482 


86.98 


(±9.76) 


0.68,0.61) 


43.1 (±7.0) 


0.794 


45.99 


(±7.81) 


0.31,0.28) 


33.38 (±5.0) 


3.867 


N/A 




0.48,0.43) 


32.51 (±5.0) 


1.169 


37.47 


(±6.28) 


0.36,0.36) 


29.11 (±5.0) 


3.053 


86.04 


(±14.85) 


0.55,0.56) 


28.3 (±6.0) 


2.375 


52.75 


(±10.6) 


1.37,1.3) 


26.01 (±8.0) 


9.32 


N/A 




0.89,0.66) 


25.52 (±6.0) 


4.122 


N/A 




1.85,1.59) 


24.75 (±9.0) 


12.036 


N/A 




1.37,1.25) 


24.66 (±7.0) 


4.008 


N/A 




1.04,0.94) 


23.09 (±7.0) 


1.31 


27.62 


(±8.02) 


0.5,0.43) 


22.97 (±5.0) 


3.046 


67.52 


(±14.43) 


1.35,1.01) 


21.34 (±7.0) 


4.063 


N/A 




0.67,0.46) 


21.15 (±5.0) 


3.922 


N/A 




0.46,0.59) 


20.72 (±5.0) 


1.093 


23.45 


(±5.46) 


1.48,1.04) 


20.17 (±7.0) 


7.618 


N/A 




1.04,1.21) 


20.12 (±7.0) 


9.547 


N/A 




1.26,1.17) 


18.88 (±7.0) 


7.405 


N/A 




1.51,1.02) 


18.67 (±7.0) 


2.846 


47.13 


(±17.25) 


1.19,1.33) 


17.94 (±7.0) 


11.011 


N/A 




0.86,1.1) 


16.43 (±6.0) 


2.461 


32.19 


(±11.18) 


0.98,0.82) 


15.95 (±5.0) 


4.225 


N/A 




0.97,1.04) 


15.01 (±5.0) 


2.891 


39.16 


(±13.9) 


1.31,1.23) 


15.0 (±6.0) 


10.808 


N/A 




0.87,1.02) 


13.8 (±5.0) 


1.722 


18.91 


(±7.06) 


0.8,1.24) 


12.99 (±5.0) 


10.729 


N/A 




0.68,0.75) 


12.95 (±4.0) 


12.372 


N/A 




0.79,0.54) 


12.94 (±4.0) 


1.073 


14.58 


(±4.88) 


0.79,0.73) 


12.11 (±4.0) 


11.09 


N/A 




1.07,1.26) 


11.96 (±5.0) 


14.397 


N/A 




1.21,0.86) 


11.93 (±5.0) 


11.554 


N/A 




0.96,1.01) 


11.09 (±5.0) 


6.779 


N/A 




0.54,0.76) 


10.71 (±4.0) 


12.977 


N/A 




0.91,0.64) 


10.3 (±4.0) 


9.316 


N/A 




0.8,0.79) 


10.08 (±4.0) 


2.345 


18.46 


(±7.65) 


0.72,0.56) 


9.36 (±4.0) 


7.982 


N/A 




0.68,0.67) 


9.27 (±4.0) 


6.277 


N/A 





Table 1. Properties of the 41 radio sources detected in the EVLA map centred on the WHDF and presented in Figure [2] Sources 
marked with a circle (o) in the left-hand column are within the region of the main lobe of the EVLA primary beam where the 
direction dependent gain is greater than 0.2. These are corrected with a primary beam model, thus the flux density values in bold 
face type represent intrinsic values, whereas the ones in plain type are apparent values. Please refer to the text for further details. 



4C 00±02) in the central region of the WHDF corresponds 
to a spiral galaxy at z = 0.046 which is also associated with 
a sub-mm LABOCA detection and an X-ray source. 
WHDF-EVLA-3: This source is a low redshift 
(z = 0.2609) galaxy. 

WHDF-EVLA-4: Near infrared imaging of this source 
shows an edge-on spiral which is clumpy in ultraviolet imag- 
ing. 

WHDF-EVLA-6: This radio source corresponds to one of 
the LABOCA-detected obscured quasars at z = 1.33 and a 
faint infrared counterpart can be seen at the centre of the 
major radio peak. The radio emission is extended and re- 
solved into two components. There is no clear optical or in- 
frared counterpart for the lower radio peak suggesting that 
the quasar may host a radio jet. Paradoxically, the opti- 



cal colours show an ultraviolet excess (UVX) while the R-H 
colours are quite red. At HST resolution the source may 
just be resolved. The redder near-infrared colours may be 
explained by host domination in these bands. Assuming all 
the R band flux comes from the QSO and R-K « 2 (e.g. Mad- 
dox et al., 2012), the QSO would have K « 20.6 implying 
K ~ 18.1 for the host galaxy. Such galaxy magnitudes are 
not uncommon at this redshift in surveys such as the K20 
survey (Cimatti et al., 2002), and as evidenced further by 
the K-band luminosity function derived from the UKIDSS 
Ultra Deep Survey (Cirasuolo et al., 2010). Note also the 
general tendency for AGN at all epochs to lie in the most 
massive galaxies (Dunlop et al., 2003). 

WHDF-EVLA-8: The second LABOCA-detected ob- 
scured quasar at z = 2.12. In the HST I band image the 
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Radio ID R-band RA R-band Dec Offset R U-B B-R R-I R-Z R-H H-K Redshift 

(J2000) (J2000) (") mag Colours within 3" diameter apertures 



WHDF-EVLA- 1 


Oh 


22m 


31. 


.86s 


+0d 


21m 


27.5s 





.8 


16 


.07 


0.04 


1 


.12 


0. 


.73 


0.89 


2. 


,19 





.32 


0.046 


WHDF-EVLA-2 


Oh 


22m 


15. 


87s 


+0d 


21m 


28.8s 





.5 


19 


.20 


0.52 


2 


.70 


0. 


.92 


1.11 


3. 


.04 


- 




0.259 


WHDF-EVLA-4 


Oh 


22m 


31. 


04s 


+0d 


22m 


24.7s 





.5 


23 


.29 


0.32 


2 


.23 


1. 


.64 


1.82 


1. 


.62 


1. 


.22 


0.961 


WHDF-EVLA-5 


Oh 


22m 


29. 


,12s 


+0d 


18m 


19.2s 





.7 


21 


.62 


-0.12 


2 


.01 


0. 


.97 


1.32 


3 


.59 





.79 


0.433 


WHDF-EVLA-6 


Oh 


22m 


24. 


83s 


+0d 


20m 


11.8s 


1 


.8 


22 


.63 


-0.72 


1 


.04 


1. 


.06 


1.33 


3 


.71 


0. 


.98 


1.32 


WHDF-EVLA-7 


Oh 


22m 


39. 


94s 


+0d 


18m 


38.6s 





.6 


22 


.44 


1.20 


3 


.42 


1. 


.67 


2.05 


4. 


.27 


0. 


.90 


- 


WHDF-EVLA-8 


Oh 


22m 


22. 


.81s 


+0d 


20m 


14.6s 





.9 


24. 


.02 


-1.25 





.79 


0. 


.71 


1.12 


3. 


41 


0. 


.85 


2.11 


WHDF-EVLA-9 


Oh 


22m 


32. 


56s 


+0d 


23m 


27.6s 


1 


.1 


19 


.31 


-0.29 


1 


.62 


0. 


.78 


1.04 


2. 


72 





.76 


0.382 


WHDF-EVLA-10 


Oh 


22m 


20. 


48s 


+0d 


19m 


0.3s 


1 


.4 


20 


.15 


-0.18 


1 


.76 


0. 


.86 


1.04 


2. 


.82 


0. 


78 


0.432 


WHDF-EVLA-12 


Oh 


22m 


11. 


08s 


+0d 


19m 


50.9s 


2 


.0 


24. 


.18 


-0.33 





.48 


0. 


.47 


- 


3. 


.85 


- 




- 


WHDF-EVLA- 14 


Oh 


22m 


24. 


99s 


+0d 


16m 


47.8s 


2 


.5 


21 


.74 


-0.26 


1 


.58 


0. 


.43 


0.57 


1. 


.56 


- 




- 


WHDF-EVLA-16 


Oh 


22m 


14. 


86s 


+0d 


21m 


4.2s 


1 


.1 


20 


.50 


0.06 


1 


.92 


0. 


.88 


2.06 


3 


16 


- 




- 


WHDF-EVLA-17 


Oh 


22m 


42. 


84s 


+0d 


19m 


54.8s 


2 


.3 


21 


.60 


0.02 


2 


.24 


0. 


96 


1.16 


2. 


92 


0. 


.75 




WHDF-EVLA-18 


Oh 


22m 


37. 


,79s 


+0d 


23m 


38.1s 


1 


.3 


17 


.39 


-0.14 


1 


.31 


0. 


.70 


0.82 


2. 


.42 


0. 


.56 




WHDF-EVLA-19 


Oh 


22m 


22. 


62s 


+0d 


20m 


43.0s 





.6 


22 


.12 


-0.27 





.87 


0. 


.56 


0.89 


2. 


.72 


1. 


.11 




WHDF-EVLA-20 


Oh 


22m 


53 


22s 


+0d 


16m 


59.8s 


4 


.3 


18 


.43 


0.06 


1 


.69 


0. 


.78 


0.80 


2. 


.57 








WHDF-EVLA-23 


Oh 


22m 


32. 


09s 


+0d 


18m 


11.1s 


1 


.3 


24 


.91 


0.16 


1 


.29 


1. 


.36 




4. 


79 





.73 




WHDF-EVLA-25* 


Oh 


22m 


36. 


,78s 


+0d 


20m 


56.5s 


1 


.5 


19 


.82 


















1. 


.01 




WHDF-EVLA-27 


Oh 


22m 


16. 


,22s 


+0d 


19m 


45.7s 


1 


.4 


21 


.39 


-0.21 


2 


.01 


0. 


.78 


1.06 


2. 


.51 








WHDF-EVLA-29 


Oh 


22m 


21. 


80s 


+0d 


21m 


47.3s 


3 


.2 


22 


.29 


-0.68 





.83 


0. 


.91 


0.66 


2. 


52 





.88 




WHDF-EVLA-32 


Oh 


22m 


28. 


24s 


+0d 


21m 


53.5s 


5 


.0 


22 


.62 


-0.17 





.97 


0. 


11 


-0.11 


1. 


.51 








WHDF-EVLA-39 


Oh 


22m 


30. 


95s 


+0d 


22m 


54.5s 


1 


.8 


23 


.77 


-0.41 


1 


.17 


1. 


.00 


1.31 


3. 


.63 


1. 


.26 





* Position and magnitude are taken from the H-band image. 

Table 2. Optical and near-infrared properties of the nearest counterpart within 5" of the positions of the radio sources in Table 
[T] Magnitude, colours and positions are taken from imaging described in Metcalfe et al. (2001) and Metcalfe et al. (2006). Spectra 
are unpublished data taken with LDSS2 on the Magellan 6.5m telescope. 



Radio ID 


LABOCA 


Comment 


R 


U-B 


870/xm 


8.4GHz 


LAB-EVLA 


Redshift 




ID 




mag 




(mjy) 


(^Jy) 


Separation 




WHDF-EVLA-6 


LAB- 11 


absorbed QSO/UVX 


22.63 


-0.72 


3.4 


48 


7."7 


1.33 


WHDF-EVLA-8 


LAB-05 


absorbed QSO/UVX 


24.02 


-1.25 


1.0 


37 


3."1 


2.12 


WHDF-EVLA-15 




blank 








28 






WHDF-EVLA-19 




merger-train 


22.12 


-0.27 




23 






WHDF-EVLA-23 




stellar - QSO? 


24.91 


0.16 




17 






WHDF-EVLA-25* 




i dropout 


19.82 






32 






WHDF-EVLA-27 




merger-double 


21.39 


-0.21 




39 






WHDF-EVLA-29 




blank 








19 






WHDF-EVLA-32 


LAB-02 


blank 






4.3 


15 


5."4 




WHDF-EVLA-39 




merger-train/UVX 


23.77 


-0.41 




18 






WHDF-EVLA- 1 


LAB-06 


spiral galaxy 


16.07 


0.04 


3.9 


201 


5."5 


0.046 



* Magnitude is taken from the H-band image. 

Table 3. Summary of Tables [T] [2] and Table 1 of lBielbv et al.l ((2012) for complete sample of 10 faint EVLA sources with S < 50/xJy 
at 8.4GHz. For completeness of the sub-mm fluxes and LABOCA-EVLA separations, WHDF-EVLA-1 is also listed, although too 
bright at 8.4GHz for inclusion in the faint sample. 



object has no nucleus and looks like an edge-on spiral, but a 
nucleus is seen both in the blue and the near-infrared imag- 
ing. Again, the colours of the quasar are anomalous, having 
U-B = -1.25, i.e. very blue, whereas R-H = 3.44 i.e. very red. 
Again, the redder near-infrared colours may be explained by 
increasing host domination in these bands. On the same as- 
sumptions as for WHDF-EVLA-6 the host galaxy could have 
K ~ 19.8 and galaxies are seen out to z ~ 2 at this limit in 
the K20 survey. 

WHDF-EVLA-9: This is a low-redshift spiral galaxy 
(z = 0.382) which harbours an optical jet and has associated 
X-ray emission. 



WHDF-EVLA-19: This object shows a stellar nucleus in 
the ultraviolet, is nearly UVX and appears among a string 
of fainter objects. It is a possible quasar or merger. 

WHDF-EVLA-23: This source has a stellar appearance 
in the HST I band images. It is very red at R-H and H-K. 
The U-B colour is blue but not UVX so it could be a z > 2.2 
QSO or a modestly reddened QSO at lower redshift, perhaps 
contaminated in the infrared by the host galaxy. 

WHDF-EVLA-25: This is a drop-out at I-band and other 
optical bands. It is only detected in H imaging, where it is 
relatively bright. 
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Figure 3. Positions and IDs of the radio sources listed in Table 
[T] and the coverage of the optical and infrared frames from which 
the thumbnails in Figure l4l have been extracted. The concentric 
grey contours show the main lobe of the EVLA primary beam. 
The outer contour is the 20% gain level and the contours have 
increments of 10%. Also shown is the area covered by the extended 
WHDF observations. 



WHDF-EVLA-27: This is a possible double object, a can- 
didate merger with relatively blue colours. 
WHDF-EVLA-39: The UVX here implies probable 
quasar. Infrared and HST I-band imaging resolves this 
source into multiple components so it could be a merger. 
Otherwise, a possible strong gravitational lens? In the lat- 
ter case, the uppermost pair of images may form a partial 
Einstein ring. 

Of the 12 sources within the ~ 5' half-power point, 9 
have optical/near-infrared counterparts within 2" and 3 are 
UVX. 4 are also LABOCA sources. Similarly, in the full area, 
of the complete sample of 10 faint (< 50/iJy) sources (see 
Table [3} 6 have optical/near-infrared counterparts of which 
3 are UVX (i.e. likely quasars) including the two absorbed 
quasar LABOCA sources. The other 3 faint radio counter- 
parts are not UVX but are only slightly less blue and likely 
to be star-forming galaxies, predominantly at lower lumi- 
nosities and redshifts. Of the 10, there are 3 possible merg- 
ers, one of which contains a UVX source. The 4 faint opti- 
cally unidentified radio sources may be either dust obscured 
quasars or galaxies. One of these sources is identified only 
at H and K and it appears resolved even at ground-based 
resolution. 

The LABOCA 870/xm fluxes of the four radio sources 
with LABOCA counterparts are also given for completeness 
in Table |3] Their high ratio of sub-mm relative to radio flux 
suggests that all four sources are likely to be dust rather 
than synchrotron dominated in the sub-mm. 



4.3 Faint radio source counts at 8.4 GHz 

Determining the counts of extragalactic sources at radio 
wavelengths has been an active area of study for several 
decades; de Zotti et al. (2010) present a review of both the 
history and the state of the art. Early radio surveys provided 
a key fulcrum in Steady State versus Big Bang cosmology 
debates in the 1950s, and observations since then have re- 
vealed much about cosmology and the evolution of radio 
sources with cosmic time, and have (particularly at higher 
frequencies) proved essential for categorising extragalactic 
foregrounds for Cosmic Microwave Background experiments. 
Source counts at the faint (sjl mjy) end of the distribution 
exhibit a turn-up. This is generally explained by the increase 
in the dominance of star-forming galaxies over AGN at these 
low luminosities (e.g. Padovani et al., 2009) and it has also 
been claimed that radio-weak AGN make a significant con- 
tribution (e.g. Jarvis & Rawlings, 2004, Smolcic et al., 2009, 
Simpson et al., 2006, 2012), although the exact nature of 
this excess remains a source of debate. While most studies 
of source counts derived from radio surveys have been con- 
ducted at L-band, the turn up persists in higher frequency 
observations including those such as the X-band observa- 
tions presented in this paper. 

The recent Absolute Radiometer for Cosmology, As- 
trophysics and Diffuse Emission (ARCADE2) experiment 
(Fixscn et al., 2011) has also piqued interest in the faint 
end of radio source populations due to the measured excess 
in the sky brightness temperature at 3 GHz (Seiffert et al., 
2011). Vernstrom et al. (2011) use source count data from 
150 MHz to 8.4 GHz to predict the contribution to the sky 
temperature background from measured source populations 
and conclude that if the ARCADE2 result is correct there 
must be an additional significant population of radio galax- 
ies at fluxes fainter than those hitherto reached by radio 
surveys. This clearly motivates the need for deeper radio 
continuum surveys. 

The source counts derived from the EVLA observations 
of the WHDF do not reach the depths needed to begin to 
address the ARCADE2 result; the counts are presented here 
for completeness. It is noteworthy however that the depth 
of the observations in this paper (2.5 fJ,Jy) is approaching 
that of the deepest 8.4 GHz observations (1.49 pJy in the 
SA13 field; Fomalont et al. 2002) which required 190 hours 
of integration time with the old VLA system. The survey 
speed advantage is brought about by the huge increase in 
available bandwidth (A*y), since the noise level in a syn- 
thesis image a oc Ais~ 0,5 . Increasing the bandwidth of the 
observation does not change the traditional limiting factor 
for high frequency survey work, which is that the sky area 
covered by a single pointing is proportional to v~ 2 where 
v is the observing frequency, but it does mean that far less 
time is required to reach a certain depth per pointing. High 
frequency surveys covering significant sky areas thus become 
far more economical. Note that since the WHDF radio ob- 
servations were taken the bandwidth available for general 
users of the EVLA has increased by a factor of 8. 

The source counts from the WHDF observation are 
shown in Figure [6] along with three other previously pub- 
lished data sets. The counts are normalized to those ex- 
pected in a non-expanding Euclidean universe for ease of 
comparison to the published data. The origins of the previ- 
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Figure 4. Near infrared / ultraviolet thumbnail images centred on the EVLA radio sources within the central area of the WHDF. The 
radio sources are presented in the right hand column. Photometric wavebands and the instrument used are noted above each column, 
radio source IDs are noted at the end of each row. The base contour level on the radio cutouts is 2cr and increases in multiples of \/2- 
Each thumbnail spans 25 arcseconds. This figure continues on the next page. 
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Table 4. Source counts and bin details for the data plotted in 
Figure [6] 



Bin 


Central flux (/^Jy) 


Bin width (fiJy) 


N 


1 


15.2 


10.4 


3 


2 


30.9 


22.0 


5 


3 


62.9 


43.0 


4 


1 


128.2 


87.5 


3 


5 


260.9 


178.1 


2 



ously published points are shown on the figure, the values of 



which were collated and listed by de Zotti et al. (2010). Only 
the 17 sources from the WHDF observations which could be 
reliably corrected for beam attenuation effects were included 
in the count. Omitting this step would clearly artificially 
bias the counts towards fainter levels. Sources were counted 
in five bins between 10 and 350 ^tJy with logarithmically 
increasing widths. Error bars are simply derived from Pois- 
son statistics. The solid angle covered by the single pointing 
within the region where the beam gain is greater than a 
factor of 0.2 is 0.0119 deg 2 for the central frequency. 

For reference, also shown on Figure [6] are the source 
counts predicted by two models down to a flux limit of 
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Figure 5. Near infrared / ultraviolet thumbnail images centred on the EVLA radio sources within the extended area of the WHDF. As 
per Figure |4] photometric wavebands and the instrument used are noted above each column. Radio source IDs are noted at the end of 
each row and the radio sources themselves are shown in the right hand column. The base contour level on the radio cutouts is 2a and 
increases in multiples of \/2. Each thumbnail spans 25 arcseconds. 



1 fiJy. The small circles show the source counts generated 
by binning all ~260 million galaxies in the full 20 x 20 deg 2 
sky area of the semi-empirical galaxy simulation of Wilman 
et al. (2008). Linear extrapolation of the 4.86 and 18 GHz 
fluxes offered by the simulation was used to generate the 
source count data. The solid line shows the model of de 
Zotti et al. (2005), solely for AGN-powered sources. The 
faint count data lie above the AGN model, illustrating the 
need for another component at faint fluxes. In the Wilman et 
al. (2008) model, this is represented by star-forming galax- 
ies (67%) and also radio- weak AGN (20%), with the lat- 
ter therefore making a comparable contribution to radio- 
loud AGN (13% j3- The apparent overestimation of the total 



model counts at the faintest limits may explain the more 
«1:1 split of AGN (both radio- loud and -weak) and star- 
forming/merging galaxies seen in Table 2 compared to the 
~2:1 ratio in favour of star- forming galaxies predicted by 
the model. Further discussion of the nature of the faint ra- 
dio source population in the WHDF is presented in Sections 
[5] and H 



5 These fractions were determined by extracting a sample of 
galaxies with fluxes between 10 and 50 fijy from 25 square degrees 
of the Wilman et al. (2008) simulation. These 100444 sources were 



then binned according to their star formation and AGN types: 
20508 radio quiet AGN; 11694 FR-I; FR-H; 1069 gigahertz 
peaked spectrum; 55894 quiescent star forming galaxies; 11279 
starburst galaxies. 
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Figure 6. Euclidean-normalized differential source counts at 8.4 
GHz. Observational data are from Windhorst et al. (1993), Fo- 
malont et al. (2002) and Henkel k. Partridge (2005), collated by 
de Zotti et al. (2010), and also from this paper. Only sources 
which have been corrected for primary beam attenuation were 
used when determining the source counts for the latter. Sources 
were placed into five bins of logarithmically increasing width be- 
tween 10 and 350 /^Jy, the lower limit being equivalent to the 
completeness limit. The error bars on the counts from this paper 
are derived from Poisson statistics. Simulated source counts are 
from de Zotti et al. (2005) and from Wilman et al. (2008). 



5 DISCUSSION 

The EVLA survey of the WHDF confirms the associa- 
tions between LABOCA sub-mm sources and X-ray ab- 
sorbed AGN. Two of the three X-ray sources associated 
with LABOCA sources are hard X-ray sources and likely 
to be cold gas absorbed, with the third X-ray source be- 
ing associated with the nucleus of a nearby galaxy and too 
faint to measure its hardness ratio. All three sub-mm+X- 
ray sources are associated with radio sources. Although the 
numbers are small, the X-ray fraction is high and so is the 
gas-absorbed fraction. Previously, Page et al. (2004), who 
studied samples of absorbed (density of Hydrogen atoms 
Nh > 10 22 cm -3 ) and unabsorbed quasars, suggested a sim- 
ilar differe nce in the sub-mm pro perties of these two classes 
of sources. iHill fc Shanks! <|2011al ) themselves found that the 
LAB OCA ECDFS sub-mm Survey (LESS, IWardlow et ail 
l201ll ) showed clear evidence for ~ 20% of sub-mm sources 
being associated with X-ray absorbed AGN. 

Models whereby absorbed AGN are used to fit the X- 



ray background with the absorbed photon s being re-rad iated 
in the infrared have b een discusssed by iGunnl (| 19991) and 



increases with hydrogen column density, the re-radiated far- 
infrared flux increases in proportion and they therefore pre- 
dict that there will be more sub-mm radiation from X-ray 
absorbed AGN. Such explanations therefore have implica- 
tions for the unified model (e.g. Antonucci, 1993) which 
holds that the differing X-ray and optical properties of ob- 
scured and unobscured AGN can be explained in terms of 
the inclination of the system with respect to the line of sight 
of the observer. The amount of sub-mm emission, originating 
in the dusty torus surrounding the system, should then be 
independent of its orientation. The differing sub-mm proper- 
ties between two classes of AGN thus cannot be explained by 
invoking the unified model. Previously it was suggested that 
sub-mm-bright AGN are being observed during a different 
evolutionary phase, whereby the AGN experiences a period 
of growth within a dusty, starforming galaxy environment 
(Page et al., 2004). It is speculated here that X-ray absorbed 
AGN might be expected to be strong sub-mm sources if the 
presence of cold, neutral gas in an AGN powered source also 
indicated the presence of cold dust. In this case the dust 
torus would need to extend up to ~1 k pc to maintain a dust 
temp erature as low as 30-35 K (e.g. iKuraszkiewicz et al.l 
l2003h . 

It is also noted that the optical and near-infrared 
colours of the two X-ray absorbed quasars are anomalous, 
being quite UVX in U-B yet quite red in R-H. If this 
is not due to variability then it might suggest that the 
optical colours are active nucleus dominated whereas the 
near-infrared colours are host dominated. Clearly this could 
favour a picture where the nuclear sight line is cold gas ab- 
sorbed as evidenced by both the hard X-rays and the optical 
narrow lines. But the dust would then have to form a more 
clumpy or toroidal structure around the nucleus to leave the 
nuclear sight-line unobscured by dust. The fact that narrow, 
high-ionisation lines like C IV are seen suggest the gas ab- 
sorption is on smaller scales than the narrow line region, 
however there appears to be little dust on the nuclear sight- 
line. This behaviour may be more characteristic of a unified 
model but such a model does not explain the basic sub-mm 
- neutral gas absorption correlation. A non-unified model 
which does explain this correlation would then also need 
to invoke a mechanism, maybe a jet, to destroy cold dust 
just along the nuclear sightline while leaving the neutral gas 
more or less in place. It is interesting that the sub-mm ab- 
sorbed X-ray AGN, WHDF-EVLA-6, does show evidence of 
a radio jet. 

The two absorbed QSOs also have radio and 
sub-mm fluxes that put them close to the FIR- 
radio correlation (e.g. Jarvis et al. 2010 and ref- 
erences therein). For example, WHDF-EVLA-08 at 



= 2.12 has S 



870/im — 
23 



6.9 

-1 



X 



10 25 WHz 



-i 



and 



iGunn fc Shanks! dT999l ) and more recently bv lHill fc Shanks! 
( 2011al ). In these models, as the absorbed X-ray fraction 



z 

Ls.4G/fz = 8.6 x 10 23 WHz -1 . An approximate con- 
version assuming dust temperature components of 35- 
60 K in the wavelength range 8-1000 /im gives 
L FIR Rilx 10 39 W = 3 x 10 12 L Q . Similarly, an approxi- 
mate conversion of the radio flux assuming a v~ x spectrum 
gives hiAGHz ~ 5 x 10 24 WHz -1 . This gives qm « 1.75 (see 
Jarvis et al 2010) compared to the average for star-forming 
galaxies of q/n « 2.2, indicating some degree of contamina- 
tion by non-thermal radio emission. 

IHill fc Shanks! (|2011al ) have argued that in a non-unified 
model the absorbed AGN may make up a significant frac- 
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tion of the bright sub-mm sources. Also, Martfnez-Sansigre 
et al. (2005) have inferred from radio observations that most 
black hole accretion is obscured, implying that the contri- 
bution of obscured AGN to the sub-mm background may be 
very significant. 

This non-unified picture where absorbed AGN are pref- 
erentially sub-mm loud is als o supported by se veral more 
recent results. For example, I Page et al.l ([2012) find that 
most of the AGN detected in the Chandra Deep Field 
North at 250 /im by Herschel SPIRE show strong X-ray ab- 
sorption, although these authors emphasised more the lack 
of far -infrared emission fro m the intrinsically bright X-ray 
AGN. ISimpson et all (|2012f l have also reported that two high 
redshift 2SLAQ quasars detected in Herschel ATLAS as 
strong far-infrared sources appear to show heavily absorbed 
Lyman-a lines. Rovilos et al. (2012, priv. comm.) using Her- 
schel PACS/SPIRE have also found that the cold dust com- 
ponent in ECDFS AGN is as significant a component in 
AGN spectra l energy distr i bution s (SEDs) as their hot dust 
components. iHickox et al.l (|2012h have also found that the 
galaxy group environment and clustering of sub-mm sources 
are complete ly consistent wi t h thos e of QSOs. Finally, a crit- 
icism of the iHill fc Shanks! i|2011al ) model is that the dust 
mass required will have to be around 5 x 10 6 Mq for L* 
QSOs, implying a cold gas mass of perhaps ~5 x 10 8 Mq. 
It is interesting to note that iMolinari et al.l (|201lT ) have de- 
tected a «100 pc ring with ~5 x 10 7 M© of gas and dust 
around the Sgr A black hole candidate in the centre of the 
Milky Way. If this was associated with a «10 43 ergs -1 AGN 
outburst then it would at least be of the right size and mass 
to produce a cold sub-mm component in the Milky Way. 

A more typical non-unified model invokes an evolu- 
tionary relation between absorbed and unabsorbed QSOs 
(e.g. Sanders et al., 1988, Fabian 1999). Here a merger trig- 
gers the birth of an obscured QSO, accompanied by a burst 
of star-formation. Once the spheroid is formed, expulsion 
of cool gas and dust quenches the star-formation, turns 
off the sub-mm emission and leaves the QSO unobscured. 
This model implies that more absorbed QSOs will be more 
FIR/sub-mm bright. Such a model would thus also be con- 
sistent with X-ray absorbed QSOs being preferentially sub- 
mm bright as obs erved in this paper. Th e main difference 
with the model of IHill fc Shanks! (|2011al ) is that there the 
dust is AGN heated wherea s in thi s other case it is star- 
burst heated. IHill fc Shanks! (l2011al ) argued that starburst 
heating leaves the similarity of the sub-mm sources' lumi- 
nosities to that of QSOs looking like a coincidence. High res- 
olution ALMA observations could also measure the size of 
the sub-mm emitting region and help differentiate between 
these possibilities. 

Herschel ATLAS results from Bonfield et al. (2011) find 
a correlation between optical QSO and FIR luminosity at 
fixed redshift as well as with redshift. This is consistent 
with the model of Hill & Shanks (2011a) since at fixed col- 
umn and redshift, higher optical/X-ray luminosity implies 
higher FIR luminosity. The amount of sub-mm emission 
from QSOs with relatively low absorption of log(Nn) ~ 20.5 
or A v w 0.1-0.2 mag, implies that 10-20% of the QSO optical 
luminosity will be re-radiated by dust, implying no disagree- 
ment with the substantial FIR luminosities for these objects 
as observed by Bonfield et al. (2011; see also Hatziminaoglou 
et al., 2010, Serjeant et al, 2010). 



For QSO dust components at higher temperatures mea- 
sured at shorter wavelengths, there may be more evidence 
for a unified model from finding correlations between AGN 
radio and IR/FIR dust emission if the radio flux is assumed 
to be orientation independent. Shi et al. (2005) found some 
indication of a correlation between 70/xm and radio fluxes 
for a sample of AGN observed with Spitzer MIPS. Previ- 
ously, Polletta et al. (2000) found that in a sample of 22 
AGN, the amount of IR/FIR dust emission was reasonably 
independent of whether the AGN were radio loud or not, 
again broadly consistent with a unified picture. Whether 
these results present a problem for either of the non-unified 
scenarios discussed above awaits further data. 

The identification of optical/near- infrared counterparts 
for 10 faint (<50 /tJy) radio sources show a mixture of ~ 
30% lower redshift, blue, starforrning/merging radio galaxies 
(3/10) and « 30% higher redshift quasars/AGN (3/10), in- 
cluding the absorbed quasars identified as sub-mm sources. 
The remaining 4/10 faint sources are optically blank, al- 
though one is detected in the near-infrared as a resolved 
source. Otherwise, there are few sources identified with red, 
early-type galaxies. Thus it might be concluded that this 
mix of sources at 8.4 GHz is not dissimilar to the mix of 
sources found i n the s ub-mm samples at 350 GHz. Indeed, 
IHill fc Shanks! (|2011bh suggested that this mix persists to 
the Herschel SPIRE bands at 350 /im, i.e. 850 GHz, with 
the AGN component then becoming increasingly less domi- 
nant at higher frequencies. The results presented above seem 
to support the idea that a strong AGN component is seen at 
all frequencies from the radio to the far-infrared, with the 
AGN component synchrotron-dominated at low frequencies 
and increasingly cold dust-dominated at higher frequencies. 



6 CONCLUSIONS 

Radio continuum observations can now achieve extreme 
depths with relatively short integration times due to the 
vastly increased correlator bandwidths that are now avail- 
able. This makes surveys much more economical, thus 
prospects are good for future, deep high-frequency surveys 
covering significant sky areas with instruments such as the 
EVLA and MeerKAT (e.g. Jarvis, 2011; Heywood et al., 
2011). There is a caveat here however: both deeper sensi- 
tivity limits and broad bandwidths conspire to increase the 
calibration complexity. Direction dependent effects (which 
are often also highly frequency-dependent) need to be prop- 
erly accounted for during calibration in order to remove the 
artefacts associated with off-axis sources. As observations 
are pushed deeper these effects become more apparent (see 
also ISmirnovll201lT ) ■ and if the science targets are extremely 
faint they are in danger of being swamped by residual cali- 
bration errors. It is fair to describe the location of the con- 
fusing source which blighted the WHDF observations as the 
'worst-case scenario' (see Appendix for details) yet the suc- 
cess with which it was removed is encouraging if future deep 
radio continuum surveys are to routinely produce images 
which are limited by thermal noise as opposed to artefacts 
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brought about by deficiencies in either the model of the sky 
or the instrument 

Applying these techniques to EVLA data, a deep 8.4 
GHz radio image of an area covering the Extended WHDF 
has been generated, and a catalogue of 41 radio sources with 
flux densities exceeding 4<r has been derived. The central, 
deepest area of the WHDF is well matched to the main lobe 
of the EVLA primary beam, and within this area 17 sources 
are detected that can have their apparent fluxes corrected 
for the beam attenuation. 

Two of these radio sources (WHDF-EVLA-6 and 
WHDF-EVLA-8) confirm the association of two high red- 
shift (z = 1.33 and 2.12 respectively) X- ray absorbed quasar s 
with the sub-mm sources detected by iBielbv et all (|2012h . 
Such sources warrant further investigation due to the uncer- 
tainty surrounding the contributions that AGN make to the 
sub-mm background. Both the sub-mm absorbed quasars 
show unabsorbed nuclear colours in the blue that would re- 
quire a mechanism to remove cold dust but not neutral gas 
from the quasar sightline. It is interesting that one of the 
sub-mm sources, WHDF-EVLA-6, appears to harbour a ra- 
dio jet. Certainly the EVLA + LA BOCA results from the 
WHDF support the previous result of lfflll fc Shanks! (|2011aT ) 
from ECDFS/LESS that X-ray absorbed AGN are signifi- 
cant sub-mm sources. This in turn may support their non- 
unified absorbed AGN X-ray background model where AGN 
dust emission is proportional to their X-ray column which 
predicts that up to ~40% of the sub-mm background may 
be due to AGN. 

The beam-corrected source fluxes are used to determine 
differential source counts which are in good agreement with 
previously published values covering the distribution down 
to «50-100 ^tJy and the semi-empirical extragalactic simu- 
lation of Wilman et al. (2008). The counts show a signifi- 
cant excess over the purely AGN-powered count model of 
de Zotti et al. (2005). The optical identifications suggest 
the WHDF faint source counts are composed of 30% high 
redshift AGN, 30% low redshift star-forming galaxies and 
40/30% optical/near-infrared blank fields, likely either to 
be dusty star-forming galaxies or the more heavily absorbed 
AGN needed to explain the X-ray background. Various in- 
dividually interesting counterparts within this faint radio 
source population have been noted. 

Future EVLA observations over a wider area of the 
WHDF will test how representative the results for the 
present area are. Deeper optical and near-infrared obser- 
vations are needed to determine further the nature of the 
population of faint radio sources responsible for both the 
turn-up in faint source counts and the excess radio sky tem- 
perature detected by the ARCADE2 experiment. 
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APPENDIX A: INTERPRETATION OF THE 
DIFFERENTIAL GAIN SOLUTIONS 

Self-calibration of a radio interferometric data set involves 
the generation of a model visibility set from assumed models 
of the sky and the instrument. Model visibilities are usually 
predicted by evaluating the radio interferometer measure- 
ment equation (RIME; e.g. Smirnov, 2011a). In the most 
commonly employed method of calibration the RIME has a 
solvable complex gain term per feed, per antenna, and a nu- 
merical algorithm is used to minimise the difference between 
the observed and the model visibility data. The best-fitting 
complex gain terms (which are functions of both time and 
frequency) are applied to correct the observed data and re- 
move the instrumental gain drifts. These corrected data are 
then imaged for deconvolution or to refine the sky model for 
further calibration. 

In the case where there is a single dominating source 
in the field (e.g. a typical observation of a phase calibrator) 
both the observed and model visibility function, and the 
associated complex gain corrections, will be dominated by 
the contribution of this source. If the dominating source is 
away from the pointing centre (as is the case for the observa- 
tions presented in this paper) then the time- and frequency- 
dependent behaviour of the complex gain solutions derived 
from self-calibration will be dominated by any instrumental 
or atmospheric effects which are specific to the direction to- 
wards that source. In the case where the target field contains 
many strong sources, each of which is subject to a direction- 
dependent effect (DDE; e.g. primary beam or ionospheric 
effects) the traditional approach of solving for a single com- 
plex gain term per antenna is insufficient if high fidelity or 
high dynamic range imaging is required. 

Continuum imaging in the presence of sources which 
are blighted by strong DDEs results in maps whose dy- 
namic range is limited by calibration artefacts rather than 
thermal noise, as is the ideal outcome. Self-calibration of 
the field will generate solutions dominated by the strongest 
source(s). Deriving a single complex gain correction is in- 
sufficient, and deficiencies in either the instrumental or sky 
model will manifest themselves as corrupted versions of the 
point spread function centred on the fainter sources, which 
deconvolution is unable to remove. 
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Figure Al. Mean amplitudes of the differential gain solutions for three different Measurement Sets (rows A, B and C). Each column 
corresponds to an antenna in the array. 



The usual approach to mitigating these effects is to em- 
ploy some form of peeling algorithm (e.g. Noordam, 2004). 
This is an iterative process whereby the sources are treated 
in order of decreasing brightness. Self-calibration is per- 
formed on a per-source basis on a visibility set which is 
phase-rotated to the position of the source in question. The 
best fitting model is computed and subtracted from the orig- 
inal data. These residuals are then used as the starting point 
for the second brightest problem source, and the process re- 
peats until an acceptable map is achieved. Aside from the 
fact that this process is somewhat unwieldy and prone to 
user error, it may have trouble converging in the scenario 
where there are multiple confusing sources of similar bright- 
ness. 

A more flexible and generic alternative now exists in the 
form of the differential gains algorithm (Smirnov, 2011b). 
This avoids the generation of intermediate data products 
by forming simultaneous solutions for the complex receiver 
gains against an all-inclusive sky model on short timescales, 
as well as solving for additional gain terms on longer 
timescales for a subset of dominating sources. The solution 
interval for the second component should be matched to the 
maximum interval over which one expects the DDE to be 
roughly constant. This maximises the signal to noise in the 
measurement, minimises the degrees of freedom in the fit, 
and effectively time-smears out contributions from the other 
sources in the field. For source subtraction purposes (as im- 
plemented in this paper) the best fitting visibility models 
for each of the sources for which differential gain solutions 
are computed are subtracted from the data and the residual 
data are imaged. 

Al Gain solutions 

The complex gain solutions themselves encode much infor- 
mation (Smirnov, 2011c). For the calibration process to be 
successful there is an implicit (and justifiable) assumption 
that (direction-dependent) gains vary smoothly with fre- 
quency and time. As mentioned in Section ^. 1 l examining the 
solutions can provide a valuable diagnostic by highlighting 
deviations from smooth behaviour. The mean amplitudes of 
the differential gain solutions for 4C 00+02 for three differ- 
ent Measurement Sets (dubbed A, B and C) can be see in 
Figure IA1I Each column corresponds to a single EVLA an- 
tenna. The plots show the mean amplitude of each solution 
as a function of time, averaged across the band. Rows A 
and B are derived from shorter 1.5 hour SBs and row C is 
derived from a 3.5 hour observation. 

The smoothness of these solutions in time is immedi- 



ately apparent. As severe as the presence of this confusing 
source was in terms of obtaining a scientifically useful map, 
the DDEs that it was subject to can be tracked by the solver, 
and despite the factor of ~100 by which 4C 00+02 was at- 
tenuated by the primary beam response there is still ample 
signal to noise to form the solutions. 

Discontinuities in the solutions (e.g. 27-B) are indicative 
of amplitude spikes on antennas that were missed during 
initial flagging. As the amplitude jumps suddenly the gain 
correction drops to compensate. 

Can the origin of the observed temporal gain variations 
be deduced? The most significant causes of DDEs in inter- 
ferometry data are the ionosphere and the primary beam 
response of the elements that make up the array. The iono- 
sphere manifests itself as a dynamic phase screen over the 
array (Interna et al., 2009) and is a significant problem for 
observing with long baseline arrays at low frequencies, how- 
ever the X-band observations presented in this paper are 
unlikely to be significantly affected by it. 

A2 The EVLA primary beam 

As mentioned in Section 12.11 it is speculated that the gain 
drifts are caused by the a combination of the structure in 
the primary beam pattern and its apparent rotation on the 
sky as the observation progresses. Figure IA2I shows a simu- 
lated complex beam pattern for the EVLA at 8.4 GHz gen- 
erated using the Cassbeam package (Walter Brisken, private 
communication). Cassbeam takes a parametrized model of 
the Cassegrain antenna optics (Brisken, 2003) and uses ray 
tracing to compute a 2x2 Jones matrix (Smirnov, 2011a 
and references therein) describing the effect of the antenna 
on the incoming radiation as a function of direction. This 
complex- valued matrix is visualised in Figure lA"2l 

The real components of the diagonal terms (LL and RR) 
show clearly the main lobe of the beam and the azimuthal 
structure beyond this induced by the antenna optics. The 
off-diagonal terms (LR and RL) effectively show the instru- 
mental polarization of the EVLA, the so-called 'beam squint' 
introduced by the fact that the two receptors sensitive to or- 
thogonal polarization modes are not co-spatial. Polarization 
will not be discussed further. 

Such beam models can also be used by the A-projection 
algorithm (Bhatnagar et al., 2008) both to predict model 
visibilities during calibration, and to apply a correction for 
known beam effects during imaging. 

It has already been noted that the confusing source in 
the observations of the WHDF was situated close to the 
first sidelobe. Rotation of the beam on the sky will therefore 
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Figure A2. Complex antenna beam patterns for all four polar- 
izations products as generated by the Cassbeam software. The 
polarization product and whether the beam pattern corresponds 
to the real or imaginary component is indicated above each frame. 
These are rendered as FITS images and fed through a MeqTrees 
simulation module to investigate the effects of primary beam pat- 
terns on various source configurations. For the diagonal terms in 
the Jones matrix (LL and RR) the base contour level is ±0.005, 
with adjacent positive and negative contours increasing in mul- 
tiples of \/2. The off-diagonal terms (RL and LR) have a base 
contour level of ±0.0003 and the same increments as the diagonal 
terms. The images span a sky area of approximately 54 arcmin- 
utes. 



cause the gain to drift according to the azimuthal asymme- 
tries. This clearly does not tell the whole story as the gains 
in Figure lAll do not exhibit similar behaviour from antenna 
to antenna. The hypothesised explanation for this is point- 
ing error: radio telescopes have varying degrees of pointing 
accuracy which has the effect of shifting the beam patterns 
shown in Figure lA"2l away from the nominal pointing centre. 
The chance positioning of the confusing source in a region 
of the beam that has azimuthal structure and steep radial 
gradients causes the gain to exhibit strong behaviour in time 
and frequency. The shifting of the beam pattern on the sky 
due to pointing error causes these variations to differ from 
antenna to antenna. The next step is to attempt a simulation 
which replicates the observations. 



A3 Simulated observation 

The CASA sm tool was used to generate a Measurement 
Set with pointing direction and frequency range consistent 
with the real observations but with a six-hour track length, 
and antenna positions matching those of the the EVLA D- 
configuration with all 27 dishes. As per the averaged real 
data 8 x 32 MHz channels were used. 

The Measurement Set was then filled with simulated 
visibilities using the Siamese framework within MeqTrees. 
The key part of this framework is the BeamSims module 
which is able to read the simulated beams presented above 
as gridded FITS images. The module performs interpola- 
tion of the beam patterns as well as sky rotation. Crucially, 
pointing errors can also be applied. Each antenna is assigned 
a random pointing offset in two orthogonal directions on the 
sky, with a value between and 10 arcseconds, a conserva- 
tive estimate of the true pointing accuracy of the EVLA 
(Vivek Dhawan, private communication). 

Two scenarios are now simulated. The first is with a 
source that is 6 arcminutes from the pointing centre, consis- 
tent with the location of 4C 00±02 in the WHDF observa- 
tions. The second is a simulation of a source at a radius of 9 
arcminutes, which places it in the approximate centre of the 
first sidelobe. As the visibilities are computed the complex 
values of the applied beam gain for each of these sources, 
per integration time, per channel, per antenna are exported 
from the RIME tree. Two channels corresponding to 8.364 
and 8.556 GHz are selected and plotted as a function of time 
in Figure IA3I The upper two rows show the amplitude val- 
ues and the lower two rows show the phases. The thickest 
line represents the lowest frequency channel. 

This simulation at least qualitatively reproduces the ob- 
served behaviour. For the source with a separation of 6 ar- 
cminutes there is large vertical scatter in the beam gains be- 
tween antennas. There are periods where gains rise as others 
fall and the variation is strongly chromatic. By contrast the 
source at a radial separation of 9 arcminutes, although it 
exhibits temporal variation due to the rotation of the beam 
pattern, exhibits much less variation between antennas. The 
source at 6 arcminutes is in a part of the beam which has 
strong gain gradients, and the pointing error exacerbates the 
disparities between antennas. 
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Figure A3. Simulated values of the complex beam gain (blue = amplitude, normalized to the maximal gain close to the beam centre; 
green = phase, in radians) as a function of time for each antenna in a six-hour EVLA simulation. Note that the phases in the r = 9' 
case have been rotated by it radians to make the plot clearer. A random pointing error with a maximum value of 10" is applied to each 
antenna, equivalent to shifting the beam patterns shown in Figure IA2I by a random amount in two orthogonal directions on the sky. For 
each complex component the plots show the behaviour for a source at a radius of 6 arcminutes from the phase centre (as per the EVLA 
observations of the WHDF) and at a radius of 9 arcminutes from the phase centre (corresponding to the approximate centre of the first 
sidelobe) . The two lines in each panel correspond to two different channels in the simulation in order to show the strong chromatic effect 
in the r = 6' case. The thick line is the 8.364 GHz and the thin line is 8.556 GHz, corresponding to the channels that span most of the 
bandwidth of the actual observations. The salient point that this plot makes is that a combination of pointing error, azimuthal structure 
in the beam and the chance positioning of a strong source in a part of the beam exhibiting strong radial gradients can result in beam 
gains with significantly different temporal and spectral behaviour from antenna to antenna. 
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